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Refined constraints on chameleon theories are calculated for atom-interferometry experiments, 
using a numerical approach consisting in solving for a four-region model the static and spherically 
symmetric Klein-Gordon equation for the chameleon field. By modeling not only the test mass and 
the vacuum chamber but also its walls and the exterior environment, the method allows to probe 
new effects on the scalar field profile and the induced acceleration of atoms. In the case of a weakly 
perturbing test mass, the effect of the wall is to enhance the field profile and to lower the acceleration 
inside the chamber by up to one order of magnitude. In the thin-shell regime, results are found 
to be in good agreement with the analytical estimations, when measurements are realized in the 
immediate vicinity of the test mass. Close to the vacuum chamber wall, the acceleration becomes 
negative and potentially measurable. This prediction could be used to discriminate between fifth- 
force effects and systematic experimental uncertainties, by doing the experiment at several key 
positions inside the vacuum chamber. For the chameleon potential V(4>) = A 4+ “ /<j> a and a coupling 
function A((f) = exp ((j>/M), one finds M > 7 x 10 16 GeV, independently of the power-law index. 

For = A 4 (l + A/</>), one finds M > 10 14 GeV. A sensitivity of a ~ 10 _11 m/s 2 would probe the 
model up to the Planck scale. Finally, a proposal for a second experimental set-up, in a vacuum room, 
is presented. In this case, Planckian values of M could be probed provided that a ~ 10 -10 m/s 2 , 
a limit reachable by future experiments. Our method can easily be extended to constrain other 
models with a screening mechanism, such as symmetron, dilaton and f(R) theories. 

PACS numbers: 


I. INTRODUCTION 


The accelerated cosmic expansion, highlighted by several 
observations such as type-la supernovae [1], the temper¬ 
ature fluctuations of the cosmic microwave background 
[2] and the distribution of large scale structures, has been 
the subject of intense research. The simplest explanation 
invoking a cosmological constant, in fair agreement with 
current observations, suffers from some fine-tuning and 
coincidence issues. A wide range of alternative models 
from various frameworks have been proposed (see e.g. [3] 
for a review), a lot of them introducing a dynamical scalar 
field coupled to matter. In order to leave observational 
signatures on the structure formation or the expansion 
history, while passing tight constraints coming from lab¬ 
oratory experiments HI], tests in the solar system (see 
[7] for a review, also [EM]) and on galactic scales mm, 
a screening mechanism appears to be fruitful by suppress¬ 
ing the fifth-force induced by the coupled scalar field in 
local environments. Chameleons mm are typical mod¬ 
els where the scalar field is suppressed in a dense medium, 
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like in the Solar system, while acting as a cosmological 
constant on sparse environment like in the cosmos at late 
time. 

Many experimental tests of the chameleons have been 
proposed so far, involving e.g. ultra-cold neutrons [151123] 
or cold atoms interferometry [23] (25] (see also BMIEQ- 
1341 ). still leaving a part of the parameter space uncon¬ 
strained. Very recently, new experiments based on atom 
interferometry, involving a source mass inside a vacuum 
chamber, have been proposed to test chameleon models 
with a high sensitivity 1351 [351 , individual atoms being 
sufficiently small to let the scalar field unscreened even 
if the nucleus is dense. The experimental setup consists 
in measuring or constraining the additional acceleration 
on individual atoms, due to the scalar field gradient in¬ 
duced by the presence of a source mass at the center of 
the chamber. Forecasts were provided in [35j and first 
experimental results of [36] claimed to rule out most of 
the chameleon parameter space. 

Nevertheless those results rely on analytical assumptions 
and have not been validated when accounting for the en¬ 
tire environment surrounding the experiment. In this 
paper, we consider a four-region model, including the 
wall of the chamber and the exterior environment, and 
we solve numerically the Klein-Gordon (KG) equation in 
the static and spherically symmetric case. The only as¬ 
sumptions made here are (1) that the scalar field reaches 
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its attractor (the field value at the minimum of its effec¬ 
tive potential) outside the vacuum chamber, i.e. in the 
air, at spatial infinity, and (2) that we have a non-singular 
C 1 solution everywhere. Our method reveals that in the 
regime where the test mass perturbs weakly the scalar 
field profile, its amplitude is controlled by the attractor 
outside the chamber instead of being related to the cham¬ 
ber size. In addition, our analysis presents the advantage 
of including the important effects of the vacuum cham¬ 
ber wall that have been neglected so far. We show that 
the scalar field profile and the resulting acceleration can 
differ by up to one order of magnitude in the weakly per¬ 
turbing regime, compared to previous analysis [351136j . 
Thus the effects of the wall cannot be neglected in a pre¬ 
cise investigation of the chameleon parameter space. In 
the thin-shell regime, analytical approximations are vali¬ 
dated to a good accuracy. Nevertheless our method high¬ 
lights order of ten percents deviations. It also allows to 
determine quantitatively the negative acceleration close 
to the walls of the vacuum chamber. 

We also propose a new experimental set-up, where the 
atom interferometer is placed inside a vacuum room. Re¬ 
fined constraints on the chameleon parameters (essen¬ 
tially the coupling function) are provided for both ex¬ 
perimental set-ups and the experimental requirements to 
exclude the chameleon parameter space up to the Planck 
scale in the future are evaluated. 

The numerical analysis has been performed for two typ¬ 
ical chameleon potentials, with varying power-laws. The 
first one (referred as Chameleon-1), V(<f) = A 4+ “/<)> a {37] 
is able to reproduce the cosmic acceleration and to fit 
supernovae data, however tests of general relativity in 
the solar system would exclude the corresponding pa¬ 
rameter space [38] ■ It is nevertheless interesting to con¬ 
sider this model as an illustrative example of field con¬ 
figurations inside the vacuum chamber that are weakly 
perturbed by the presence of the central source mass. 
Laboratory experiments also probe other regions of its 
parameter space, even if they are not of direct inter¬ 
est for cosmology. The second considered potential (re¬ 
ferred as Chameleon-2) has an additional cosmological 

constant [55] . V{(f>) = A 4 exp ^A a /(fi a ^J — A 4 (1 + A a /<j) a ) 

where we assume A = A for keeping only one additional 
parameter, a being fixed. For both models we take an 
exponential coupling function A(<f) = exp(</>/M). 


The paper is organized as follow. In Sec. [TT] we briefly 
remind the equations of motion and introduce the con¬ 
sidered models. The experimental setup is described in 
Sec. [HT] and the numerical strategy is detailed in Sec. m 
Analytical results are reminded in Sec. |V| a nd are com¬ 
pared to the numerical results in Sec. |VI| for the two 
considered models, refined constraints being established 
on their parameters. We finally discuss our results and 


draw some conclusions and perspectives in Sec. VII 


II. THE MODELS 


We start from a generic action for modified gravity mod¬ 
els involving a non-minimally coupled dynamical scalar 
field 0, written in the Einstein frame, 



S->) 2 -^) 
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with R , the scalar curvature, k = 8i r/m^, m p \ being 
the Planck mass, the matter fields, V(</>) a general 
potential and a general coupling function. 

Varying this action with respect to <\> gives the KG equa¬ 
tion 




dV din A 

d(f) d <f> 
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where T is the trace of the energy-momentum ten¬ 
sor Tp,, = -(2/y/^g) (dSm/dg^). We introduce 
the stress-energy tensor for a perfect fluid in the Jor¬ 
dan frame in order to consider conserved quantities, i.e. 

= 0, the tilde denoting Jordan frame quantities. 
Energy density p and pressure p in both frames are then 
related by gS] 

P = A 4 {cj))p , (3) 

P = A i {(j))p , (4) 


so that in the weak field regime the KG equation in the 
static and spherically symmetric case becomes 



dV e ff 

d(j) 


dV 
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dA 
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(5) 


in which we have introduced an effective potential V e s, 
the prime denoting a radial coordinate derivative. We 
work here in the non-relativistic limit with negligible met¬ 
ric potentials and pressure. In Table[IJ the two considered 
inverse power-law chameleon potentials and the coupling 
function A{(j>) are specified. Equations of motion for both 
models are identical since only the derivative of the po¬ 
tential dV/d(f) contributes to the KG equation. However 
viable A values are different. For the Chameleon-1 model, 
A must be fixed by the supernovae best-fit and obeys 

to m, 

190 — 47 

logA[GeV] ~ —j— , (6) 

4 + a 


in order to reproduce the acceleration of the Universe ex¬ 
pansion even in the presence of the non-minimal coupling 
[381 . In the case of Chameleon-2, A is the cosmological 
constant value A ~ 2.4 meV while a is an independent 
variable. 


III. EXPERIMENTAL SETUP 

Laboratory experiments measuring the acceleration in¬ 
duced by a test mass can be used to probe and constrain 
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Potential V(rf>) 

Coupling function A((j >) 

Parameters 

Minimum of V e g (</ min ) 

Mass sq. at minimum (m^, in ) 

Chameleon 1 

A 4 +“/</>“ 


(a, A), M 



1/(Q + 1) 

a(l + «)A 4 +“(^|^)^ 

Chameleon 2 

a 4 (i+a %n 


a, A, M 



!/(“+!) 

a(l + a)A 4 +“(^|^)^ 


TABLE I: Characterization of the effective potential for the models of interest: potential V(4>), coupling functions A(<j>), the 
model parameters, the minimum of the effective potential and the corresponding mass. Dependent parameters appear into 
brackets. Values in the last two columns are valid as long as A(</) ~ 1 (an assumption no longer valid for some models, see 
discussion in Sec VI l. 


modifications of gravity. As a reminder, the chameleon 
field is screened in high density environments while it me¬ 
diates long-range force in sparse ones. Therefore atomic 
particles in a ultra-high vacuum chamber can mimic cos¬ 
mos conditions. In the first experiments, the test mass 
was located outside the vacuum chamber [2U [35], an ex¬ 
perimental setup which is not ideal given that the cham¬ 
ber wall screens the fifth force on the atoms. New experi¬ 
ments have been proposed in [35., SB] where the test mass 
is located inside the vacuum chamber, which improves 
the constraints on the acceleration due to the scalar held. 
Here we focus on a recently proposed atom interferome¬ 
try experiment [36] where one takes advantage of matter- 
wave properties of cesiunr-133 atoms in a Fabry-Perot 
cavity. When an atom absorbs/emits a photon, it re¬ 
coils with a momentum p = hk, with k the wavenumber 
of the absorbed/emitted photon. So, one can reproduce 
the equivalent of a Max-Zelmder interferometer for cold 
atoms with three light pulses using counter-propagating 
laser beams. Atoms are initially prepared in a hyperhne 
state F = 3 and stored in a 2 dimensional magneto¬ 
optical trap. A first light pulse splits the matter-wave 
packet in two hyperhne state F = 3 and F = 4 and 
gives an impulse of hk e g to the atoms. The effective 
wavenumber k e g depends on the two counterpropagating 
beam wavenumbers. The probability of hyperhne tran¬ 
sition can be controlled by the intensity and duration of 
both laser beams. The second pulse reverses the rela¬ 
tive motion of the beams like the mirror of Max-Zehnder 
interferometer and the third pulse acts like a beam split¬ 
ter which allows overlap of partial matter wave packets. 
Because of the recoil of the atoms, the phase difference 
between the two arms of the interferometer A</> is a func¬ 
tion of the acceleration a of atoms, 

A(/> = fceffaT 2 , (7) 

where T ~ 10 ms in general, is the time interval be¬ 
tween two pulses. To alleviate some systematics effects, 
counterpropagating laser beams are reversed and the alu¬ 
minum sphere can be positioned in two places: a near 
and a far positions (the test mass surface is respectively 
located 8.8 mm and 3 cm far from the atoms), which al¬ 
lows to disentangle the contribution from chameleon force 
to Earth’s gravity. One measurement consists thus of 
four interference fringes, corresponding to reversed coun¬ 
terpropagating laser beams and both positions of the test 
mass. Using this setup, the acceleration induced by the 


Ra 

Radius of the test mass 

lcm/5.1 x 10 13 GeV _1 

L 

Radius of the chamber 

10cm/5.1 x 10 14 GeV _1 

Rw 

Wall thickness 

lcm/5.1 x 10 13 GeV _1 

m A 

Test mass 

11.3g /6.7 x 10 24 GeV 

PA 

Test mass density 

1.2 x 10 _17 GeV 4 

Pw 

Wall density 

3.5 x 10 _17 GeV 4 

Pv 

Vacuum density 

5.0 x 10 _35 GeV 4 

Patm 

Air density (P a tm) 

5.2 x 10“ 21 GeV 4 


TABLE II: Fiducial experimental parameters, corresponding 
to the setup of [36]. 

chameleon has been excluded up to 

a exp < 5.5 /im/s 2 at 95% C.L.. (8) 

The experimental setup proposed in [ 35] is similar, except 
that they plan to use cooled rubidium atoms launched 
in a small fountain located 1 cm far from the test mass. 
Our numerical simulations can be easily adapted for such 
a configuration. 

Details of the considered experimental setup are reported 
in Table |ll] The size and density of the central mass, 
the geometry of the chamber and the vacuum density 
are those of [35, 35] ■ In addition we consider the thick¬ 
ness and density of the vacuum chamber walls, as well as 
the exterior density. In Figjl] we draw the experimen¬ 
tal setup considered in our numerical simulations. The 
four regions are labeled by their densities 1 : (1) the test 
mass made of aluminum (pa), (2) the vacuum where the 
acceleration due to the chameleon is measured (p v ), (3) 
the wall of the chamber (p w ) made of stainless steel, (4) 
the exterior of the chamber, mostly filled by air at atmo¬ 
spheric pressure (p a tm)- 


IV. NUMERICAL STRATEGY 

Analytical approaches have been considered so far [35, 
136] , which are valid under some assumptions like negligi¬ 
ble chamber wall effects. Therefore, numerical methods 


1 In the following of the paper, p refers to the density in the Jordan 
frame. 
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FIG. 1: Outline of the atom-interferometry experiment, sim¬ 
ulated by a four-region model including the source mass, the 
vacuum chamber, its walls and the exterior environment. In 
light gray, the near and far positions where the acceleration 
on atoms is measured (note that we consider a fixed source 
mass to keep spherical symmetry whereas in the real experi¬ 
mental setup the source mass is moved [36]). 


are useful to validate and refine analytical results, by in¬ 
cluding the effects due to the experimental setup, like the 
thickness and the density of the wall as well as the exte¬ 
rior environment. In the future, numerical results will be 
also helpful to study more realistic situations where the 
vacuum chamber is not exactly spherical or cylindrical. 
We consider two methods for solving the KG equation 
©: a singular and multipoint boundary value prob¬ 
lem (bvp) solver with unknown parameter and a non¬ 
linear bvp solver implementing up to sixth order a mono- 
implicit Runge-Kutta method with an adaptative mesh 
refinement, working in quad precision 2 . In the latter case, 
the density in the four regions was made continuous by 
considering arctan profiles with negligible widths. 

We take the minimal assumption which states that the 
scalar field is settled to its attractor at spatial infinity, i.e. 
0oo = ^min(patm) as [56]. Then, the asymptotic scalar 
field profile is obtained by linearizing the KG equation 
up to first order around spatial infinity, 

ft' H— 4> — A4 2 (<^ — 0oo) i (9) 

r 

with A4 2 = d 2 I4ff/d^ 2 |^ ^ , which admits the Yukawa 
profile solution (A4 2 > 0) 

Ce~ Mr 

<t> = 4>oo H- , (10) 

r 


2 For this purpose we have used the Matlab function bvp4c which 
deals with singular bvp’s and a modified version of the mirkdc 
bvp solver with adaptative mesh in Fortran. 


with C the constant of integration. Since the KG equa¬ 
tion is of second order and the parameter C is to be de¬ 
termined, three boundary conditions are needed. They 
are provided by the regularity condition on the scalar 
field derivative at the origin ft(r = 0) =0 and by the 
asymptotic behavior of cj) and ft given by Eq.(10) at the 
end of the integration interval. For the multipoint bvp 
method, the continuity of (j> and ft are imposed at the 
interfaces of each region (6 conditions) while the profile 
is guaranteed to be continuous for arctan profiles of den¬ 
sity. The density and size of each region are reported in 
Table HU The two numerical methods have been checked 
to be in agreement with each other. Their applicabil¬ 
ity to the various regimes and their limitations in the 
deep thin-shell regime will be discussed in Sec. [Vi] We 
already point out that this numerical method enables 
to properly account for the effect of neighboring matter 
on the chameleon fields and can be easily generalized to 
other experiments, possibly more sensitive (in the limit 
of spherical symmetry). 

V. ANALYTICAL APPROACH 

In this section we reproduce the main steps of [55] and de¬ 
rive analytically the chameleon held profile in the spheri¬ 
cally symmetric and static regime for a two-region model 
(the source mass and the vacuum chamber). In the 
next section, the validity of the various assumptions will 
be analyzed and the analytical approximations will be 
compared to the exact numerical results, for the two 
chameleon potentials of Table [I] For the sake of simplic¬ 
ity, we assume in this section that a = 1. 

Assuming A(ft) = 1, the held value at the minimum of 
the effective potential written in Eq. ([b]), and the held 
mass around it, are respectively given by 



The case where the effect of A((f>) is important will be 
discussed in Sec. m For a two-region model the density 
p is either the source mass density pa or the density in 
the vacuum chamber p v . 

Four different regimes can be identihed, depending on 
whether the held reaches the effective potential minimum 
or not: (1) the held does not reach the minimum of the 
effective potential in any region, (2) the held reaches the 
minimum in the vacuum chamber but not in the source 
mass, (3) the held reaches the minimum in the source 
mass but not in the vacuum chamber, (4) the held reaches 
the minimum both inside the test mass and the vacuum 
chamber. Cases (1) and (2) were referred as the weakly 
perturbing regime in Ref. [35] . whereas (3) and (4) were 
referred as strongly perturbing. Below we consider those 
four cases separately, as in Ref. (36]. In principle, one 
should also distinguish between the cases where the held 
reaches 0 m i n inside the chamber wall, or not. When low¬ 
ering M, depending on the central mass density and size, 
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on the chamber wall density and thickness, </> m in can be 
reached first inside the central mass or inside the chamber 
walls. Nevertheless, for the considered experimental set¬ 
up, the wall and the central mass have similar densities 
and sizes, and so those two cases will not be distinguished 
in the following. 


1. <j)(r = 0) ^ 0min(pA) and (/>(Ra < r < L) -/= 0 min (pv) 


Within the test mass the field does not reach the attrac¬ 
tor that is the minimum of the effective potential. Since 
p v < Patm < Pa > the second term in the effective poten¬ 
tial dominates, T4ff — (f>p A /M. The KG equation can be 
solved inside the mass imposing that the field profile is 
regular at the origin, which gives 


(f> = D + 


m\r 2 
8ttMR\ ’ 


( 12 ) 


where D is an integration constant that can be fixed by 
matching <fi and q{/ to the field solution in the vacuum 
chamber at r = Ra- Inside the vacuum chamber the 
field does not reach the attractor value. Let us denote 
cj) b g the value that would take the field at the center of 
the chamber in the absence of the source. Then one can 
consider an harmonic expansion of the potential 


9 

Td 

V e f[(4>) CS v e g((/) hg ) H-^(0 - 4> bg) 2 , (13) 

higher order terms being subdominant. One can solve 
the KG equation assuming that the field profile decays 
at infinity. This gives 

4 > { r ) = </>b g + “e- m ^ r • (14) 

r 

Note that at r = Ra , one has mg g RA 1 for typical 
experimental parameters and thus 4>{Ra) — 4>b g + cz/Ra- 
After matching, one finds the held profile in case (1) 


<p W (r) = </> bg 


m A 

8itRaM 



Ra c ~ m ^ r ) 
r ) 


0(r - Ra) 


r) 

(15) 


Color 

M [GeV] 

a <t>!g {near) 

0-4-19 if 

ar) 

Chameleon-1, weakly perturbing: Figs. 2 

3 

Blue 

Green 

Red 

Light blue 

10 13 

10 15 

10 17 

10 19 

1.3 x 10 1 
1.3 x 1CT 3 
1.3 x 10 -7 
1.3 x 10“ n 

2.8 x 10° 
2.8 x 10 -4 
2.8 x 10 -8 
2.8 x 10" 12 

Chameleon-1, thin-shell: Figs. 
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Blue 

Green 

Red 

Light blue 

10 s 

10 9 

10 10 

10 11 

5.8 x 10 9 
5.2 x 10 s 

1.9 x 10 7 
2.5 x 10 s 

L4 x 7 10 s 
5.7 x 10 6 
—4.4 x 10 6 
5.5 x 10 4 

Chameleon-2, thin-shell : Figs. 

9 

10 


Blue 

Green 

Red 

Light blue 
Purple 
Beige 

10 14 

10 15 

10 16 

10 17 

10 18 

10 19 

5.2 x 1CT 7 
5.2 x 10 -8 
5.2 x 10 -9 

5.2 x 10~ 10 

5.3 x KT 11 
4.6 x 1(T 12 

1.5 x lCT 8 
1.5 x 1(T 9 
1.5 x 10 _1 ° 
1.5 x 10" 11 
2.4 x 10" 12 
6.8 x 10" 14 


TABLE III: Properties of the numerical scalar field and ac¬ 
celeration profiles for the two models in the different regimes. 


the scalar held reaches the minimum of the effective po¬ 
tential inside the vacuum chamber wall. This assumption 
is actually not valid in case (1) because p w ~ pa, and be¬ 
cause the wall thickness is about the radius of the test 
mass. So in most of the parameter space corresponding 
to case (1), the scalar held does not reach its attractor 
inside the wall. As a result, 4>b g is better approximated 
by </>min(Patm)- Numerical results will highlight the ef¬ 
fects of the chamber wall on the scalar held profile. Even 
if the background held value has no effect on the accel¬ 
eration itself, this result is important because it changes 
the region in the parameter space in which case (1) ap¬ 
plies: it is extended to lower values of M, as developed 
thereafter. 

The analytical held profile and the induced acceleration 
ciif, = d r (j)/M have been plotted on Figs. [2] and [^respec¬ 
tively for various values of M reported in Table |Ilf| 

The acceleration induced by the scalar held gradient in¬ 
side the vacuum chamber is well approximated by 


where 0 is the Heaviside function. Therefore the effect 
of the mass is to deepen the held prohle, by a quantity 
3 toa/(87t.RaM) < <j> b g at r = 0. By dehnition, case (1) 
is valid as long as \(j>b g — ^^(r = 0)| <C </>b g . Outside the 
mass, the difference \4>g g — <p\ decreases like oc 1/r for re¬ 
alistic experimental conhgurations where the exponential 
decay factor can be neglected. 

A subtlety arises in the evaluation of <pb g , which in 
Ref. [35] was either the attractor in the vacuum, either 
related to the chamber size 3 , under the assumption that 


p v is much lower than the wall density p w where the field was 


“* “ 4 0 + mi «) ' (16) 

Since m\, g r -C 1 for realistic laboratory experiments, the 
acceleration is independent of A and thus one can con¬ 
strain directly the value of M. This is the reason why, 


assumed to reach its attractor 0 m i n (p w ) - Thus the first term of 
V eff in Eq. ([13]) dominates the KG equation inside the chamber, 
which can be solved to get </>bg as a function of the size of the 
vacuum chamber. However, behind this calculation is hidden the 
assumption that the field reaches </> m i n (pw) in the wall, which is 
not valid in case (1) in most of the parameter space. 
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FIG. 2: Numerical and analytical scalar field profiles (respec¬ 
tively solid and dashed lines) in absolute value for various M 
listed in Tabic |Th| The numerical profiles obtained when low¬ 
ering the wall density to p w = 5 x 10~ 19 < p\, are drawn 
in dotted line, in order to illustrate that the wall density can 
perturb more or less importantly the field profile. Vertical 
lines mark out the four regions: the test mass, the chamber, 
the wall and outside the chamber. 



1 o~ 2 IQ' 1 10° 10 1 


Radial coordinate r(L) 

FIG. 3: Numerical and analytical profiles (respectively solid 
and dashed lines) of the acceleration a^/g with g the Earth 
gravitational acceleration, for M values listed in Table |III| 
Vertical lines mark out the four regions (test mass, chamber, 
wall and exterior). 


as we will show in the following, the power-law of the 
potential has no effect on the acceleration as long as 
\A{(j))- 1|<1. 


a. (f) min (p A ) and 0 bg = 0 min (pv) 

When the size of vacuum chamber is larger than the char¬ 
acteristic distance over which the field reaches the mini¬ 
mum of the potential, that is when 


L > 


1 

tttmin (Pv) 


/ A 5 M 3 \ 1/4 

V 4 pi J 


(17) 


the field profile is still governed by Eq. ©. However 
the value of 0b g is now simply 0 m i n (Pv)- In the case 
of the bare chameleon potential V((j> ) = A 5 /0, one has 
A ~ 2.6 x 1CT 6 GeV in order to reproduce the late- 
tirne accelerated expansion of the Universe. For typi¬ 
cal vacuum densities and chamber sizes, e.g. those re¬ 
ported in Table |TTJ one finds that this regime would 
occur when M < 1CV 6 GeV. This does not corre¬ 
spond anymore to the weakly perturbing regime requiring 
0b g > m A /(4TrR A M), which gives M > 2 x 10 9 GeV in 
our fiducial experimental setup. In the case of the poten¬ 
tial V(0) = A 4 (l + A/0), A ~ 10 -12 GeV is the cosmo¬ 
logical constant. It results that the field in the chamber is 
expected to reach 0 m i„(p v ) only if M < 10 5 GeV. There 
again this is far from the regime where the test mass per¬ 
turbs only weakly the field, valid when M > 10 2 ° GeV, 
i.e. in the super-Planckian regime. 


s. m — 0 min (Pa) and 0(i?A < r < L) ^ 0 min (p v ) 


In case (3) the field reaches 0a = 0min (pa) inside the 
test mass. One can define a radius S such that 0(5) = 
0 a( 1 + e) with 0 < e <C 1. For S < r < R\, the density 
term dominates in V e g and the solution of the linearized 
KG equation is given by 


= D 


C 

r 


m\r 

8nMR 3 A 


(18) 


which is the same as Eq. (12) but with a non-vanishing 


integration constant C. Outside the test mass, the field 
still obeys to Eq. (14). After matching 0 and 0' at r = 


Ra and 0 at r = 5, the integration constants a, D and 
C can be fixed. The resulting field profile in case (3) 
reads E5] 
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0 ( 3 ) ( r ) _ < 


0A , 



0 A + 

m A / 3 

3 S 2 r + 2 S 3 

8tt R\Mr V 

0bg - 

m A c -m bx r / 

vi)’ 

4-7T Mr \ 


r<S , 

S < r < R\, 
r > Rai 


(19) 


with the radius 


S = R A \ 1- 


87rMi?A^b g 


3toa 


( 20 ) 


being such that one has typically (i?A — S)/Ra 1, 
corresponding to the thin shell regime. The induced ac¬ 
celeration is well approximated (mb g l?A <C 1) by 


_ m A 
^ 4:TrM 2 r 2 


1 - 


Ri 


Rj\ dbg 
Mr 2 


( 21 ) 


and contrary to case (1), it is related to the value of 
0b g . If the wall is sufficiently large, then the field reaches 
4>min(Pw) and so the calculation of 0b g in Ref. [35] is valid, 
giving 


0 bg ^ 0.69 (A 5 L 2 ) 


1/3 


( 22 ) 


for a spherical chamber. Compared to the case (1), the 
induced acceleration therefore does not depend only on 
M but also on A and on the size of the vacuum chamber 
L. When A is set to the cosmological constant and L 
to the fiducial value reported in Table [TT] one finds that 
the experiment of [36j constrains the coupling parameter 
down to M ~ 10 15 GeV. The above calculation does not 
involve the power-law index a (apart indirectly via mb g , 
but there is no effect in the limit mb g r <C 1). Therefore 
it is expected that the predictions are independent of a, 
as long as | A(<j>) — 1| <C 1. 

Analytical field profiles and induced accelerations for case 
(3) are represented on Figs. E] and [5] for the bare potential 
V(4>) = A 5 /(f> (Chameleon-1), and on Figs. [9] and 10 for 
the potential V(0) = A 4 (l + A/0), for several values of 
M reported in Table [TTT| Those are found to be in good 
agreement with the numerical results, except close to the 
wall where important deviations are found. 

In the strongly perturbig regime, the reliability of the 
theory is questionable. Indeed the quantum corrections, 
either in the matter and the chameleon sector must re¬ 
main small. Most of the parameter space reachable by 
the experiment proposed in ([36] belongs to this regime. 
Following ([35] the underlying instabilities are harmless 
and the classical analysis is trustable, keeping in mind 
that quantum corrections can become large at very small 
scales. However since the aim of this paper consists of 
modelling how the environement can affect the analyti¬ 
cal results derived for the classical field, we also provide 
numerical forecasts in the questionable strongly perturb¬ 
ing regime. Nevertheless we did not explore the deeply 


strong regime but focus on the transition between the 
two regimes, where the numerical computations allow to 
follow the smooth evolution of the field and acceleration 
profiles whereas analytical assumptions break. Our com¬ 
putations show that the analytical estimations are recov¬ 
ered once in the strong regime, and that they are quite 
reliable, at least classically. The underlying quantum as¬ 
pects are beyond the scope of this paper. 


4 ■ 0(9) — 0min (Pa) QTld 0 bg — 0min(Pv) 


In case (4) the field profile is governed by Eq. (19) since 
the field reaches the effective potential minimum at the 
center of the test mass. However, as long as the condi¬ 
tion Eq. (17) is satisfied, 0bg = 0 m in (p v )- For the bare 
potential V(0) = A 5 /0, case (4) takes place when M < 
10~ 3 GeV, whereas for the potential V(0) = A 4 (l-(-A/0) 
one needs M < 10 GeV in order to be in the strongly per¬ 
turbing regime inside the test mass. Therefore the case 
(4) is irrelevant for values of A motivated by cosmology 
and realistic experimental configurations. 


VI. FOUR-REGION MODEL: NUMERICAL 
RESULTS 

A. Weakly perturbing regime 

We discuss the weak field regime for Chameleon-1 ex¬ 
clusively since only the case (3) of the previous section 
(strongly perturbing regime) is relevant for values of M 
below the Planck scale in the Chanreleon-2 model. Al¬ 
lowing superplanckian values, one actually would recover 
the regime where the field is only weakly perturbed by 
the central mass, but the induced acceleration would be 
far too low for being observable with future experiments. 
On Fig. [2} |0 — 0oo| is represented for various M values of 
the Chameleon-1 model, a = 1 and A = 2.6 x 10~ 6 GeV 
being fixed, which corresponds to case (1) discussed in 
the previous section. Inside the test mass the scalar field 
is constant but the numerical profile roughly differs by 
two orders of magnitude compared to the analytical ap¬ 
proximation. This difference is induced by the effect of 
the wall, which enhances |0 — 0oo| in the vacuum cham¬ 
ber: the wall tends to stabilize the scalar field and gives 
it a kick right outside the wall shell. Note that the im¬ 
portance of the effect depends on the wall density and 















FIG. 4: Numerical results (solid lines) and analytical ap¬ 
proximation (dashed lines) for the scalar field profile of the 
Chameleon-1 model, in the strongly pert-urbing (thin shell) 
regime, for A = 2.6 x 10 -c GeV and values of the coupling 
M listed in Table [Til] Differences between the two-region and 
four-region models are non-negligible inside the chamber, es¬ 
pecially at the vicinity of the wall. Vertical lines mark out 
the four regions (test mass, chamber, wall and exterior). 



1CT 1 10° 

Radial coordinate r(L) 


thickness. By setting p w <C pa (or by reducing the wall 
thickness), one tends to recover the analytical profile. 
Outside the chamber, the scalar field follows the Yukawa 
profile as imposed by the asymptotic behavior in Eq. (101. 
The acceleration a^jg = d r (j>/{Mg) with g the Earth 
gravitational acceleration, is plotted on Fig. [3| The gen¬ 
eral behavior of the acceleration profile obtained numer¬ 
ically does not differ significantly from the analytical 
approximation. However the chamber wall affects the 
amplitude of the profile with a difference growing up 
to one order of magnitude at the vicinity of the wall. 
This result illustrates how important it is to take into 
account the four regions modeling the experiment, in 
the weakly perturbing regime, in the view of establish¬ 
ing accurate constraints from atom-interferometry exper¬ 
iments. Given the experimental constraint on the ac¬ 
celeration a exp /g < 5.6 x 10 -7 [36], we find that the 
Chameleon-1 model is excluded at 95% C.L. for a = 1 
and M < 7 x 10 16 GeV (see Table III). An experiment 
controlling systematics to probe a^/g < 10“ 12 would rule 
out the model up to the Planck scale. This is close to the 

as a reachable sensitiv- 


io- 


given m 


value a^/g • 
ity. 

On Fig. [6} the parameter space of ( a,M) is explored 
for the acceleration measured in the near position (see 
Sec. |III[ ). Devia tion s between near and far positions are 
tiny (see Table III for an order of magnitude). A uni¬ 
versal behavior with respect to a is observed for M > 
10 17 GeV. Given the current bounds on a^/g, we thus 
find that the coupling parameter M is constrained iden¬ 
tically for Chanreleon-1 models, independently of the a 
power-law parameter. Future experiments will not be 
able to distinguish between the various power-law poten¬ 
tials. The universal behavior is broken at lower values of 
M because the assumption | A(cf>) -— 1| <gc 1 is not valid 
and the analytical approach cannot be trusted anymore. 


B. Strongly perturbing regime 

Probing the deep thin-shell regime, i.e. when ( Ra,w — 
5a, w ) •C 5 a, w (thin-shell radius of the test mass or of 
the chamber wall), is very challenging numerically. Up 
to some point, it is nevertheless possible to track the so¬ 
lution and to check the validity of the analytical estima¬ 
tions, typically using mesh refinement methods. The nu¬ 
merical treatment also allows to probe the smooth tran¬ 
sition between the weakly and strongly perturbed cases. 


1. Chameleon-1 


FIG. 5: Numerical results (solid lines) and analytical approx¬ 
imation (dashed lines) for the acceleration \a^\/g profile, for 
the same model and parameters as in Fig. [4] The numerical 
profile for the four-region model shows that from the middle 
of the chamber to the wall, the acceleration becomes negative 
and increases in magnitude. Vertical lines mark out the four 
regions (test mass, chamber, wall and exterior). 


Even if one has already predicted analytically that the ac¬ 
celeration would be excluded, the field and acceleration 
profiles have been also computed for parameters corre¬ 
sponding to the strongly perturbed regime, referred as 
Case (3) in the Sec. [V| Those are represented on Figs. [4] 
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FIG. 6: Forecast for the normalized acceleration a^/g mea¬ 
sured in the ’’near” position, i.e. 8.8 mm far from the test 
mass, for various M and a. Dotted line represents the gravi¬ 
tational acceleration due to the test mass. 


Color 

Pa [GeV 4 ] 

Pw [GeV 4 ] 

Blue 

1.0 x 10“ 2 ° 

1.0 x 10“ 2 ° 

Green 

2.5 x 10~ 2 ° 

2.5 x 10 _2 ° 

Light blue 

5.0 x 10” 2 ° 

5.0 x 10“ 2 ° 

Purple 

7.5 x 10~ 20 

7.5 x 10“ 2 ° 

Beige 

5.0 x 10" 19 

7.5 x 10“ 2 ° 

Red 

1.2 x 10” 17 

7.5 x 10“ 2 ° 


TABLE IV: Densities inside the test mass p\ and the wall 
p w for the numerical scalar field and acceleration profiles of 
Figs. 0and[8] 


and [5] and compared to analytical predictions for several 
values of M. Our numerical method also allows to probe 
the transitory regime 

As expected given that pa 1$ Pw with similar test mass 
radius and wall thickness, when lowering M, the field 
reaches first the potential minimum </> m i n (p w ) inside the 
wall, and then (f> m i n (PA ) within the test mass, over a 
very thin radius. Inside the vacuum chamber, one ob¬ 
serves that the field roughly reaches the amplitude of 
<j) b g given by Eq. (22), which validates the calculation 
of j35j . In the vicinity of the chamber wall, however, the 
acceleration changes its sign and becomes negative, with 
a comparable magnitude with the acceleration close to 
the source mass. This effect could be helpful experimen¬ 
tally to discriminate between a signal of modified gravity 
and systematic errors, by performing measurements of 



FIG. 7: Numerical results (solid lines) and analytical ap¬ 
proximation (dashed line) for the scalar field profile of the 
Cha mele on-2 model, for various p\ and p w reported in 
M = 10 1 ' GeV being fixed. 


Tab. 


IV 



FIG. 8: Numerical results (solid lines) and analytical ap¬ 
proximation (dashed line) for the acceleration ap/g profile 
of th e C hameleon-2 model, for various pA and p w reported in 
M = 10 1 ' GeV being fixed. 


Tab. 


IV 


the acceleration at several key positions of the chamber. 
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2. Chameleon-2 


For the Chameleon-2 model and the considered experi¬ 
mental set-up, it has been impossible to track numerically 
the thin-shell regime up to the point where the acceler¬ 
ation would have been large enough to be observed in 
laboratory experiments. Nevertheless, the field and ac¬ 
celeration profiles are represented on Figs. [7] and [8j for 
M = 10 17 GeV and increasing values of p w and p A - The 
attractor field values within the test mass and the wall 
are reached progressively and the field variations at the 
borders between the four regions become more steep, as 
expected given that (Ra,w ~ S a , w )/Ra,w cx Mp~^ w R~^ 2 w 
[see Eq. m- I n the case M = 10 17 GeV, the attractor 
is reached inside the test mass for p A — 5 x 10 -2O GeV 4 , 
i.e. about 1000 times lower than the aluminum den¬ 
sity, whereas inside the wall, it is reached for p w ~ 
7.5 x 10 _2O GeV 4 . This slight difference is explained by 
the fact that the central test mass has a diameter two 
times larger than the wall thickness. 

Inside the vacuum chamber, the analytical estimation is 
roughly recovered in the first half of the chamber. Once 
in the thin-shell regime, one can also observe that the 
field and acceleration profiles inside the chamber are in¬ 
dependent of the wall and mass densities, except at their 
immediate vicinity. Therefore, in the deep thin-shell 
regime, the scalar field and acceleration both at the near 
and far positions of the interferometer do not depend on 
the test mass and wall densities and sizes, neither on the 
exterior environment. In order to obtain the numerical 
solution inside the chamber, down to low values of M , 
one can therefore use the trick to set the wall and mass 
densities high enough to be in the thin-shell regime but 
low enough for the field profile to be numerically tractable 
through the borders of the four regions. 

The field and acceleration profiles have been calculated 
numerically and compared to the analytical results, for 
several values of M and A ~ 2.4 meV. These are repre¬ 
sented on Figs. [9] and 10 As expected the profiles have 
the same behavior than for the Chameleon-1 model (see 
Figs. [4] and [5|. Close to the test mass, one recovers the 
analytical estimation but one can nevertheless notice dif¬ 
ferences higher than 20%. 

Close to the wall, the acceleration becomes negative, and 
its amplitude reaches values comparable to the acceler¬ 
ation at a position close to the test mass, which is a 
potentially measurable prediction that could be useful to 
discriminate between experimental systematic effects and 
an acceleration induced by the presence of some scalar 
field. 

In conclusion, we find that the atom-interferometry ex¬ 
periment of .36] already excludes values of the coupling 
parameter M < 10 14 GeV at 95% C.L. Moreover, if 
the experimental sensitivity could be reduced down to 
a^/g ~ 10“ 11 (as it is claimed to be feasible in [35]), the 
model would be probed nearly up to the Planck scale. 
Finally, note that the typical field values reached inside 
the chamber are too low to induce large deviations from 


A((f>) — 1, which implies that our results are roughly in¬ 
dependent of the power-law index a. 


C. Chamber geometry effects 


The numerical method used throughout this paper takes 
into account the effects of the chamber geometry, in the 
limit where the vacuum chamber is spherical. Explor¬ 
ing various chamber size and wall density, we propose 
to consider the possibility to realize the same atom in¬ 
terferometry experiment in a vacuum room in order to 
make the test of M values up to the Planck scale pos¬ 
sible in a near future. The largest vacuum rooms have 
a radius larger than R = 10 m and their walls made of 
concrete are sufficiently large for the field to reach inside 
its effective potential minimum. One can thus neglect 
the exterior of the chamber (see Sec. VIBl. The vacuum 
room can sustain a vacuum around 10 -D Torr (we assume 
p v = 5 x 10“ 31 GeV 4 ), low enough to prevent </> bg to reach 
its effective potential minimum in vacuum. 

Numerical field and acceleration profiles are reported on 


Figs. |TT] and 12 respectively. Assuming as before pa = 
1.2x 10 -17 GeV 3 , it results that a test mass of 1 cm radius 
only enables to probe the regime where the field does 
not reach <f >a inside the test mass (see dashed green lines 
on Figs. 11 and 121, the acceleration being thus poorly 


constrained. However, provided that the test mass radius 
is larger (e.g. Ra = 3.3 cm), the strongly perturbing 
regime is reached and the acceleration is large enough to 
be measurable in a near future for M of the order of m p \. 
As a result, for M = m p 1 , {a^l/g = 2.4 x 10 _1 ° at 8.8 mm 
from the surface of t he t est mass (the previously called 
near position in Sec. Ill) while \a<f,\/g = 5.7 x 10~ 10 for 
M = 0.1 rn.pi. In comparison, the test mass of 1 cm gives 
rise to |a^,|/g = 1.7 x 10 -11 for M = m p \. 

Similarly to what was obtained in Sec. |VI B the thin shell 
regime cannot be tracked numerically if the wall density 
is of the order of the concrete p ~ 10” 1 ' GeV 4 . But 


one can safely consider lower values of p w (see Fig. 11) 


without any significant change of the results inside the 
vacuum room. 


VII. CONCLUSION 

The chameleon screening mechanism is able to suppress 
the fifth-force induced by a scalar field in locally dense 
environment, while allowing the scalar field to be respon¬ 
sible for dark energy on large astrophysical scales and 
thus to affect significantly the large scale structure forma¬ 
tion. This makes chameleon theories good candidates for 
explaining dark energy and for being testable by future 
cosmology-dedicated experiments, such as Euclid g3] or 
the next generation of giant radio-telescopes dedicated to 
21 cm cosmology [44]. Chameleon theories are also well 
constrained by local tests of gravity in the solar system, in 
the galaxy, as well as in laboratory. Recently it has been 
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FIG. 9: Numerical results (solid lines) and analytical ap¬ 
proximation (dashed lines) for the scalar field profile of the 
Chameleon-2 model, in the strongly pert-urbing (thin shell) 
regime, for A = 2.4 x 10“ 12 GeV and values of the coupling 
M listed in Table |Ul| The test mass, wall and exterior den¬ 
sities have been adapted for making the profile numerically 
tractable, with no effect inside the vacuum chamber (apart 
at t he imm ediate vicinity of the borders), as explained in 
Sec. VIB 2 The ratios M/p were kept constant (with the 
same value than for the red curve of Fig. [TJ, which fixes 
the thin-shell radius, apart for M = 10 18 GeV (purple) and 
M = 10 19 GeV (beige) for which only the wall density was 
adapted. Noticeable deviations from the analytical estima¬ 
tion are observed inside the chamber, due to the wall effects. 
Vertical lines mark out the four regions (test mass, chamber, 
wall and exterior). 



FIG. 10: Numerical results (solid lines) and analytical ap¬ 
proximation (dashed lines) for the acceleration \a$\/g profile, 
for the same model and parameters as in Fig. [9] Strong dis¬ 
crepancies are observed between the four-region (numerical) 
and the two-region (analytical) models. Vertical lines mark 
out the four regions (test mass, chamber, wall and exterior). 


FIG. 11: Numerical scalar field profiles for a vacuum room 
(L = 10 m). The green dashed curve is obtained for a test 
mass of Ra =1 cm with M = m p (p w = 2.5 x 10 -21 GeV 4 ) 
while the blue and the red ones are obtained for Ra = 3.3 cm 
with M = 0.1 x m p (p w = 2.5 x 10 -22 GeV 4 ) and M = m v 
(p w = 2.5 x 10 -21 GeV 4 ) respectively. We only consider a 
three regions model, neglecting the effect of the exterior of 
the vacuum room (see discussion in Sec. VIBI. 



FIG. 12: Numerical acceleration profiles \a^/g\ for a vacuum 
room (L = 10 m). The green dashed curve is obtained for a 
test mass of Ra =1 cm with M = m p (p w = 2.5 x 10^ 21 GeV 4 ) 
while the blue and the red ones are obtained for Ra = 3.3 cm 
with M = 0.1 x m v (p w = 2.5 x 10 -22 GeV 4 ) and M — m v 
(p w = 2.5 x 10~ 21 GeV 4 ) respectively. We only consider a 
three regions model, neglecting the effect of the exterior of 
the vacuum room (see discussion in Sec. VIBI. 
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proposed to use an atom-interferometry experiment to 
constrain chameleon models with an unprecedented ac¬ 
curacy by probing the acceleration induced by the pres¬ 
ence of the scalar field on cold atoms. The experiment 
is realized inside a vacuum chamber in order to reduce 
the screening effect, and a central mass is used to source 
some field gradient. Forecasts were calculated in [33 and 
a first experimental setup was build and used to estab¬ 
lish new constraints on chameleon model [37T] . However 
the calculations of the field and acceleration profiles rely 
on several approximations, and until now did not fully 
consider the effects of the vacuum chamber wall and of 
the exterior environment. 

The purpose of this work was to validate and refine previ¬ 
ous calculations, by using a numerical approach consist¬ 
ing in solving the Klein-Gordon equation in the static and 
spherically symmetric case for a four-region model rep¬ 
resenting the central source mass, the vacuum chamber, 
its wall, and the exterior environment. Three bound¬ 
ary conditions are imposed: the field must be regular at 
origin and reach the minimum of the effective potential 
with a Yukawa profile, at large distance in the exterior 
environment. Our method allows to probe the transition 
between the regime where the central source mass only 
weakly perturbs the field configuration, and the thin-shell 
regime where the field inside the central mass and in¬ 
side the chamber walls reaches the minimum of the ef¬ 
fective potential over a very small distance. Two typical 
chameleon potentials were considered, in inverse power- 
laws and allowing varying powers, as well as a standard 
exponential form for the coupling function. 

In the weakly perturbing regime, it is found that the 
chamber wall enhances significantly the scalar field in¬ 
side the vacuum chamber and reduces the induced ac¬ 
celeration, by up to one order of magnitude compared 
to previous analytical estimations and with a maximal 
effect close to the wall. 

Going to the thin-shell regime, for our fiducial exper¬ 
imental setup, the field reaches the attractor inside the 
chamber wall and the exterior environment becomes thus 
irrelevant. However, for reasonable value of the induced 
acceleration, the field inside the vacuum chamber does 
not reach the minimum of the effective potential and is 
instead related to the size of the chamber, as first no¬ 
ticed in [35]. Our analysis refines the field and acceler¬ 
ation profiles in the chamber and highlights noticeable 
deviations from the analytical estimation, which is nev¬ 
ertheless roughly recovered close to the central test mass. 
Close to the chamber wall, the acceleration becomes neg¬ 


ative, with a magnitude similar to the one close to the 
central mass. We argue that this prediction could be 
useful to distinguish between systematic effects and fifth- 
force effects which should be maximal and opposite close 
to the central mass and to the wall , and should vanish 
roughly at the middle-distance between the test mass and 
the walls. 

Refined constraints have been derived on the coupling 
parameter M from the atom-interferometry experiment 
of [HQ. For the chameleon potential V(cf> ) = A A+a /cf> a 
and a coupling function A(</>) = exp (c/>/M), one finds 
M > 7 x 10 16 GeV, independently of the power-law. For 
the bare potential V{4>) = A 1 2 3 4 (l + A/</>), we find that 
M > 10 14 GeV. We have also confirmed that a future ex¬ 
periment reducing its sensitivity down to a ~ 10~ lo m/s 2 
would be able to rule out most of the parameter space of 
the latter model, nearly up to the Planck scale. 

Finally, we have proposed to realize a similar atom- 
interferometry experiment inside a vacuum room. The 
density inside such rooms is low enough for the field pro¬ 
file and the induced acceleration to depend only on the 
size of the room. If the room radius is larger than about 
10 meters, we find that the chameleon model could be 
probed up to the Planck scale. Nevertheless, further work 
is needed to implement realistic non-spherical geometries 
of the room (or of the vacuum chamber). 

We conclude that numerical results will be helpful in the 
future in order to establish accurate bounds on various 
modified gravity models. In particular, the effects of 
the vacuum chamber wall and its exterior environment 
cannot be neglected. Our numerical method is easily 
extendable to study other forms of the field potential 
and other modified gravity models requiring a screen¬ 
ing mechanism, such as the symmetron, dilaton and f(R) 
models. An investigation of the symmetron model is in 
progress and should be released soon. Finally it can be 
easily adapted to other experiments. 


VIII. ACKNOWLEDGMENTS 

We warmly thank Christophe Ringeval, Clare Burrage, 
Holger Muller, Justin Khoury, Benjamin Elder and 
Philipp Haslinger for useful comments and discussion. 
The work of S.C. is partially supported by the Return 
Grant program of the Belgian Science Policy (BELSPO). 
S.S. is supported by the FNRS-FRIA and A. F. is par¬ 
tially supported by the ARC convention No. 11/15-040. 


[1] A. G. Riess et al. (Supernova Search Team), Astron.J. 
116, 1009 (1998), astro-ph/9805201. 

[2] P. Ade et al. (Planck) (2015), 1502.01589. 

[3] L. Amendola and S. Tsujikawa, Dark energy: theory and 
observations, Cambridge University Press (2010). 

[4] E. Adelberger (EOT-WASH Group) (2002), hep- 


ex/0202008. 

[5] A. Upadhye, Phys.Rev. D86, 102003 (2012), 1209.0211. 

[6] D. Kapner, T. Cook, E. Adelberger, J. Gundlach, B. R. 
Heckel, et al., Phys.Rev.Lett. 98, 021101 (2007), hep- 
ph/0611184. 

[7] C. Will, EAS Publ.Ser. 30, 3 (2008). 



13 


[8] V. Anastassopoulos et al. (CAST) (2015), 1503.04561. 

[9] D. Scott and M. Rees, Mon.Not.Roy.Astron.Soc. 247 
(1990). 

[10] R. Pourhasan, N. Afshordi, R. Mann, and A. Davis, 
JCAP 1112, 005 (2011), 1109.0538. 

[11] B. Jain, V. Vikram, and J. Sakstein, Astrophys.J. 779, 
39 (2013), 1204.6044. 

[12] J. Khoury and A. Weltman, Phys.Rev.Lett. 93, 171104 
(2004), astro-ph/0309300. 

[13] J. Khoury and A. Weltman, Phys.Rev. D69, 044026 
(2004), astro-ph/0309411. 

[14] S. S. Gubser and J. Khoury, Phys.Rev. D70, 104001 
(2004), hep-ph/0405231. 

[15] P. Brax, C. van de Brack, A.-C. Davis, J. Khoury, 
and A. Weltman, Phys.Rev. D70. 123518 (2004), astro- 
ph/0408415. 

[16] D. F. Mota and D. J. Shaw, Phys.Rev.Lett. 97, 151102 
(2006), hep-ph/0606204. 

[17] D. F. Mota and D. J. Shaw, Phys.Rev. D75, 063501 
(2007), hep-ph/0608078. 

[18] T. Jenke, G. Cronenberg, J. Burgdorfer, L. Chizhova, 
P. Geltenbort, et al., Phys.Rev.Lett. 112, 151105 (2014), 
1404.4099. 

[19] H. Lemmel, P. Brax, A. Ivanov, T. Jenke, G. Pignol, 
et al., Phys.Lett. B743, 310 (2015), 1502.06023. 

[20] P. Brax and G. Pignol, Phys.Rev.Lett. 107, 111301 
(2011), 1105.3420. 

[21] P. Brax, G. Pignol, and D. Roulier, Phys.Rev. D88, 
083004 (2013), 1306.6536. 

[22] A. N. Ivanov, R. Hollwieser, T. Jenke, M. Wellenzohn, 
and H. Abele, Phys. Rev. D 87, 105013 (2013), URL 
http://link.aps.org/doi/10.1103/PhysRevD.87. 
105013 

[23] P. Brax and C. Burrage, Phys.Rev. D83, 035020 (2011), 
1010.5108. 

[24] J. Fixler, G. Foster, J. McGuirk, and M. Kasevich, Sci¬ 
ence 315, 74 (2007). 

[25] G. Lamporesi, A. Bertoldi, L. Cacciapuoti, M. Prevedelli, 
and G. Tino, Phys.Rev.Lett. 100, 050801 (2008). 

[26] A. Shih, Phys. Rev. A 9, 1507 (1974), URL http:// 
link.aps.org/doi/10.1103/PhysRevA.9.1507 

[27] A. Shih and V. A. Parsegian, Phys. Rev. A 12, 
835 (1975), URL http://link.aps.org/doi/10.1103/ 


PhysRevA.12.835 

[28] A. Anderson, S. Haroche, E. A. Hinds, W. Jhe, and 

D. Meschede, Phys. Rev. A 37, 3594 (1988), URL http: 
//link.aps.org/doi/10.1103/PhysRevA.37.3594 

[29] C. Sukenik, M. Boshier, D. Cho, V. Sandoghdar, and 

E. Hinds, Phys.Rev.Lett. 70, 560 (1993). 

[30] F. Baumgartner, R. J. Sewell, S. Eriksson, I. Llorente- 
Garcia, J. Dingjan, J. P. Cotter, and E. A. Hinds, Phys. 
Rev. Lett. 105, 243003 (2010), URL http://link.aps. 
org/doi/10.1103/PhysRevLett.105.243003 

[31] D. M. Harber, J. M. Obrecht, J. M. McGuirk, and E. A. 
Cornell, Phys. Rev. A 72, 033610 (2005), URL http: 
//link.aps.org/doi/10.1103/PhysRevA.72.033610 

[32] M. Kasevich and S. Chu, Phys. Rev. Lett. 67, 
181 (1991), URL http://link.aps.org/doi/10.1103/ 
PhysRevLett.67.181 

[33] A. D. Cronin, J. Schmiedmayer, and D. E. Pritchard, 
Rev.Mod.Phys. 81, 1051 (2009). 

[34] P. Brax, C. van de Brack, A.-C. Davis, D. F. Mota, and 
D. J. Shaw, Phys.Rev. D76, 124034 (2007), 0709.2075. 

[35] C. Burrage, E. J. Copeland, and E. Hinds, JCAP 1503, 
042 (2015), 1408.1409. 

[36] P. Hamilton, M. Jaffe, P. Haslinger, Q. Simmons, 
H. Muller, et al. (2015), 1502.03888. 

[37] B. Ratra and P. Peebles, Phys.Rev. D37, 3406 (1988). 

[38] A. Hees and A. Fuzfa, Phys.Rev. D85, 103005 (2012), 
1111.4784. 

[39] P. Brax, C. van de Brack, A.-C. Davis, J. Khoury, 
and A. Weltman, Phys.Rev. D70, 123518 (2004), astro- 
ph/0408415. 

[40] T. Damour and K. Nordtvedt, Phys.Rev. D48, 3436 
(1993). 

[41] C. Schimd, I. Tereno, J.-P. Uzan, Y. Mellier, L. van Waer- 
beke, et al., Astron.Astrophys. 463, 405 (2007), astro- 
ph/0603158. 

[42] A. Upadhye, W. Hu, and J. Khoury, Phys. Rev. Lett. 
109, 041301 (2012), 1204.3906. 

[43] L. Amendola et al. (Euclid Theory Working Group), Liv¬ 
ing Rev.Rel. 16, 6 (2013), 1206.1225. 

[44] P. Brax, S. Clesse, and A.-C. Davis, JCAP 1301, 003 
(2013), 1207.1273. 


